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ANALYTICAL SOLUTION FOR THE WIND-DRIVEN CIRCULATION 


IN A LAKE CONTAINING AN ISLAND 
by Marvin E. Goldstein and Richard T. Gedney 
Lewis Research Center 

SUMMARY 

An analysis has been carried out to determine analytically the effect of an island on 
the wind driven currents in a shallow lake (or sea). A general analysis is developed that 
can be applied to a large class of lake and island geometries and bottom topographies. 
Detailed numerical results are obtained for a circular island located eccentrically or 
concentrically in a circular lake with a logarithmic bottom topography. It is shown that 
an island can produce volume flow (vertically integrated velocities) gyres that are com- 
pletely different from those produced by a normal basin without an island. These gyres 
in the neighborhood of the island will produce different velocity patterns, which include 
the acceleration of flow near the island shore. 

INTRODUCTION 

Lake currents are very important in determining what chemical and biological 
changes will take place in a lake since they control the distributions of sediment and pol- 
lution which have been introduced. It is, therefore, important to be able to predict these 
currents. 

The main source of motion in the oceans and lakes are the winds. When the winds 
act on the surface of a deep body of water, they set up a circulation which consists of a 
top and bottom Ekman (bovtndary) layer separated by a geostrophic (inviscid) core. For 
lakes , in which the depth is much greater than the thickness of the friction or Ekman 
layers, the usual Ekman dynamics (i. e. , the bottom stress assumed proportional to be 
the geostrophic velocity) can be used. For any of the Great Lakes, the use of Ekman 
dynamics is questionable since they all have shallow shore regions of considerable 


extent. In Lake Erie, for moderate winds, the thickness of the Ekman layers are com- 
parable to the depth over much of the lake and therefore the use of Ekman dynamics is 
not valid. 

The necessary extension of the Ekman analysis to the case of a shallow lake has 
been given by Welander (ref. 1), and that theory has been used here. The solution for 
the surface displacement is first obtained. From this, the horizontal velocities are cal- 
culated as a function of depth and horizontal position in the lake. In the present analysis, 
only the motions caused by steady winds are considered. Gedney (ref. 2) has numer- 
ically calculated the velocities for Lake Erie using Welander’s shallow sea formulation. 
Comparison of these calculations with measurements in Lake Erie show that steady -state 
currents actually do occur. 

Goldstein, Braun, and Gedney (ref. 3) developed a method for obtaining product so- 
lutions to a certain class of partial differential equations. This class of equations in- 
cludes Welander 's equation for the wind-driven circulation in shallow lakes for a large 
class of bottom topographies. The method was used to reduce the problem of finding so- 
lutions to Welander’s equation for a closed body of water to the problem of solving an or- 
dinary differential equation. To illustrate the method, the complete analytical solution 
for an elliptically shaped body of water with a particular choice of bottom topography 
when the depth is greater than one-half the Ekman friction depth was given. In this anal- 
ysis the method developed in reference 3 is extended to include the case when an island 
is included. 

Previously, islands have been analyzed numerically in total studies, which included 
effects such as irregular mainland boundaries that have masked the Island effect. This 
analysis is being performed to determine just the local effect due to an island so that a 
fundamental understanding is obtained. To date no such island analyses have been per- 
formed for a shallow lake. 


BASIC EQUATIONS 

In the present analysis the basic approximations are that the water density is con- 
stant, the vertical eddy viscosity is independent of depth but dependent on wind velocity, 
the pressure is hydrostatic, and the lateral friction and nonlinear acceleration terms can 
be neglected. The horizontal scale used to jiistify neglecting horizontal friction and non- 
linear acceleration is of the order of the width of the lake. For regions close to the 
shore, usually of the order of 1 kilometer, the horizontal friction and nonlinear accelera- 
tion terms are important. Therefore, the lake must be much larger than 1 kilometer and 
the currents calculated herein are ’’large scale” ciirrents that are valid in the interior 
region away from the shore areas. 
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These assumptions reduce the momentum equations to two equations containing as 
unknowns the horizontal velocities and the surface slope. The appropriate boimdary con- 
ditions for these equations are a no -slip condition at the lake bottom and a specified 
shear stress (due to the wind) at the air -sea interface. These equations and boundary 
conditions can be solved analytically to give the velocity as a function of depth with the 
stirface wind stress and stirface slope as parameters. 

By vertically integrating the momentmn equations and by using the vertically inte- 
grated continuity equation, one can obtain a single equation for the stir face displacement: 


^xx + ^yy + 



( 1 ) 


where ^ is the surface displacement of the sea or lake, h is the depth of the lake, x 
and y are the horizontal, locally Cartesian coordinates, d is the Ekman 

friction depth, v is the coefficient of vertical eddy diffusivity, and f^ is the Coriolis 
parameter, which introduces the effect of the Earth's rotation. In addition. 
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( 5 ) 



(6) 


where g is the acceleration due to gravity and and are, respectively, the x 
and y components of the wind stress divided by the density of water. 

Equation ( 1) is a slightly rearranged form of that given by Welander (ref. 1) . For 
details of the derivation of equation (6) and for a discussion of the approximations in- 
volved in the derivation, see Gedney (ref. 2). For instance, Gedney has shown for Lake 
Erie that the approximations used, except for the approximation of a constant eddy vis- 
cosity, induce at the most a small (of the order of 10"^) or very local error in the cal- 
culations. It is shown in reference 3 that, when h > id, the coefficients in this equation 
can be approximated quite closely by much simpler functions of h/d. Thus, 
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The X and y components of the horizontal velocity U and V, respectively, can 
be found once the solution to equation (1) is known by using the relation (given in ref. 1) 


(U + iV) 


II 

■(1 -i)7rr'^l 

sinh 

(1 +i)7^| 


gd J 

cosh 

(1 + i 

dj 



cosh 7 t( 1 + i) — 
d 


0yy 


( 10 ) 


where ^ is the vertical distance below the water surface and where 


- tY + 


( 11 ) 


The X and y components of the total volume flow Qj and Qg, respectively, can 
also be fotmd from the surface displacement by using the relations 


^ = — / Ud^=J- ^ (E? - F? )] 

;d gd 27 t L d y J 


gd 

Q2^c _ ^c 


gd gd jf 




( 12 ) 


where 


sin/27T + sinh/2ir - 
d/ 


F = 


X \ u/ \ d/ ^ j 

2t 7 — C0S^27T — ^ + COSh^27T — ^ 


(13) 


When h > i d, equation (12) can be approximated closely by the simpler relation 
(ref. 3) 


27TQf^ 

gd 


,T-[a + « 


for !l>l 
d 2 


(14) 


where we have put 


Q = Qj + iQg 


(15) 
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and 



+ la 


( 2 ) 


or, using equations (9) and (11) and collecting terms, 


(16) 


ffRi i _ jl for il>i (17) 

gd J d 2 

The boundary conditions for equation (1) are obtained by specifying the normal vol- 
ume flow at the coasts. Thus, if n^ and Ug are the x and y components, respec- 
tively, of the outward drawn unit normal n to the boundary, then for any closed body of 
water the boundary conditions for equation (1) are obtained by substituting equation (12) 
(or for the case where (h/d) > (1/2) eq. (14)) into the relation 

Q^n^ + Q2n2 = 0 (18) 


Class of Depth Distributions for Which Equations Can Be Solved 

We shall now suppose that there exists a nonconstant harmonic function u of x and 
y and an arbitrary function H of u only such that the depth distribution h can be ex- 
pressed in the form 


J.H[u(x,y)] (19) 

d 27T 

This is a fairly general functional form, and it will be possible, for any one of a 
large number of lakes and seas, to choose the functions u and H in equation (19) so 
that the depth distribution is approximated closely by a relation of this type. 

When relation (19) is substituted into equations (1) and (2), we find that 

+ Cyy + 2 [q!(u)u^ - ^(u)UyJCj^ + 2[o'(u)Uy + i3(u)uj^y = f (20) 

where 
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a 


sin H sinh H ^ 

(cos H + cosh H)(sin H - sinh H) du 


^ = i sinh H + sin H dH 
2 cosh H + cos H du 


(21) 


f H 


_ cosh H + cos H 


sinh H - sin H I 9x 


. 9x 0y . 


( 22 ) 


or when h > i d the coefficients a and /3 are given to a close approximation by 
2 


a 0 


2 du J 


for h>l 


(23) 


Now, let V be the harmonic conjugate of u. Then the function 


w = u + IV 


(24) 


is an analytic function of the complex variable 

z = X + iy 

and the fimctions u and v satisfy the Cauchy -Riemann equation. Upon taking u and v 
as new independent variables for equation (20) and introducing a new dependent variable 
Z defined by 


Z 


^|fo 


sinh H - sin H 


cosh H + cos H 


or by 


Z = ? f or h > i d 

2 


(25) 


(26) 


(if the approximations (23) are used) we find, after using the Cauchy -Riemann equations 
to simplify, that equation (20) becomes 

Zuu + V + WZ^ - {[a(u)]’ + [a(u)f}z (27) 
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where 




/cosh H + cos H 

Dd 


1— 1 

y sinh H + sin H 

[°^x 

j 

Idwl 


When h > i- d the transformed equation can be approximated by 
2 



( 28 ) 


(29) 


The details of the algebraic manipulations involved in obtaining equations (27) to (29) can 
be fotind in reference 3. 

We are interested in solving equation (1) for a closed lake or sea that contains a sin- 
gle island. Thus, we shall seek solutions to equation (1) in the doubly connected region 
shown in figure 1 where the outer boundary (shore line) is “if and where the inner (is- 
land) boundary is 

We shall suppose that the bottom topography of the basin can be approximately de- 
scribed by a function of the form (19) and that the depths along the shore lines ^ and y 
are both constants . Equation (19) now shows that the harmonic function u must also be 
constant along both if and ^ . Without any loss of generality we can always redefine the 
functions H and u so that u = 0 on We shall denote the constant value of u on y 

by u . Again, without loss of generality, it is always possible to arrange matters so 
that u > 0. Now, consider the anal 3 rtic function w defined in terms of u by equa- 
tion (24) and put 


T = e^ (30) 

Then, T is an analytic function within the miiltiply connected region of figure 1 and since 

[T ( = e" 

u 

it follows that |T | = 1 for z on and |T ( = e ® = constant for z on 5^ Hence the 
mapping 


z -T 
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transforms the multiply connected region of figure 1 into the interior of the concentric 
annular region in the T-plane shown in figure 2, with the outer boundary mapping into 
the (inner) unit circle and the inner boimdary y mapping into the outer circle with radius 


u„ 


Ro = e 


( 31 ) 


On the other hand, as long as the boundaries of the regions shown in figure 1 are 
sufficiently smooth, the Reimann mapping theorem for doubly connected regions guar- 
antees that there always exists an R^ for which such a mapping exists (ref. 4). 

The mapping w = In T transforms the cut annular r^ion in the T-plane (see fig. 3) 
into the rectangular region in the w -plane shown in figure 4. Thus, u can always be 
chosen so that the anal 3 rtic function w, defined in equation (24), maps the cut region oc- 
cupied by the body of water shown in figure 5 into the rectangular regions in the w-plane 
shown in figure 4. The cut occurs along a line of constant v. 

Now we have transformed equation (1) into equation (27) or, when h >id, into 
equation (29) with u and v as the independent variables. Thus the problem of solving 
equation ( 1) for the lake configuration shown in figure 1 has been transformed into the 
problems of finding a solution to equation (27) or, alternatively, when h > — d equa- 
tion (29) in the rectangular region in the w-plane shown in figure 4. Although this sim- 
plifies the shape of the region in which the problem is to be solved, the most important 
simplification is due to the fact that all the coefficients in equations (27) and (29) are 
functions of only one of the independent variables. As will be shown subsequently, this 
will allow an analytical solution to the problem to be found. 

Before this is done, however, it is first necessary to transform the boundary condi- 
tion (18) in the physical plane into a boundary condition in the w-plane. 

To this end, notice that since u = 0 on the boundary and u = u on y the unit 

O 

normal vector n is given by 


n 



on ^ 


n = 



on y 


Hence, it follows from the Cauchy-Riemann conditions and the definition of the derivative 
of an analytic function that the components of n can be written as 


9 


and 





dw dz 
dz dw 


dw 

dz 


dw 





J 


(32) 


Thus, it follows from equations (12), (16), and (19) that boundary condition (18) can be 
written as 


Re [ct + H(E + iF)fii + i = 0 

L \3x 3y/JVdw/ 


for z on and for z on y 


After using the Cauchy -Riemann equations and the chain rule for partial deriva- 
tives, this becomes 


HF 


— - HE ^ = Re ) on u = 0 and u = u for -n v ^ n 


3v 3u 


\dw/ 


It now follows from equations (25), (6), (13), and (19) that Z must satisfy the following 
boundary conditions. 


Pq^u + + SqZ r^°^(v) for u = 0; -ir < tt 

Ps^u + ‘^s^v ®s^ " r^^Vv) for u = Ug; -tt < v < 


(33) 


where 



p _ sinh H(0) - sin H(0) 
° cosh H(0) + cos H(0) 


sinh H(u„) - sin H(u ) 
s s 

cosh H(u„) + cos H(u ) 
s s 


_ jj^Q^ _ sinh H(0) + sin H(0) 
^ cosh H(0) + cos H(0) 


% = H(Ug) 


sinh H(u ) + sin H(u ) 
s s 


cosh H(u„) + cos H(u ) 

S 0 


S„ = 


sin H(0) sinh H(0) 


° [cos H(0) + cosh H(0)]^ 
sin H(Ug) sinh H(Ug) 


[cos H(Ug) + cosh H(Ug)j 


2 Vdu 


u=u„ 


(34) 


and 




i 




sinh H(0) - sin H(0) 
cosh H(0) + cos H(0) 


- sin H(u ) 


V sinh H(Ug) - si 
cosh H(u ) + cc 


,g, , cos H(Ug) 


(^e 


Re 


a(^ 

Vdw 


u=0 


■(-) 

Vdw/ 


u=u_ 


(35) 


When h > id, equation (33) can be replaced by the approximate boundary conditions 


S ^U ^S ^ V 


(°>(V) 

for u = 0; 


for u = Ug; 


for ->i 
d 2 


(36) 


where now 
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and 


Pq = Ps = 1 


% - H(0) - 1 


for ll>i 
d 2 


qs = H(u^) - ij 




> 

y^°\v) = Re 

\dw/ 

u=0 


r riT 

> 

y(s)(v) = Re 

k-) 



\dw/ 



d 2 


( 37 ) 


(38) 


Since equation (27) is elliptic, it is still necessary to specify boundary conditions 
along the sides v = ±ir of the rectangular region shown in figure 4. These conditions 
follow from the fact that ^ and its normal derivative must be continuous across the cut 
in the physical plane (see fig. 5). Since the opposite sides of this cut map into the lines 
V = +7T and V = -n, we must require that 


and 



(39) 


Or in view of equation (26) these can be replaced by the condition 



when the approximate equation (29) is to be used. 
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We have now shown that the surfe.ce dispfecement for a sea or lake of arbitrary 
shape can be fovind by solving the differential equation (27) in the interior of the rectan- 
gular region shown in figure 4 subject to the boundary conditions (33) and the periodicity 
conditions (39). Or, if h/d > 1/2, the surface displacement is approximately obtained 
by solving equation (29) in the interior of the rectangular region in figure 4 subject to the 
boundary conditions (36) and the periodicity conditions (40). 


Reduction to Ordinary Differential Equations and Summary 


In order to solve the boundary value problem, we put 


S2j^(u) = 




Z(u, v)e“^^^ dv n = 0, ±1, ±2, 


(41) 


Then it follows from the theory of Fourier series that this transform can be inverted to 
obtain 


OO 

Z(u,v) = ^ J2jj(u)e^*^'" 

n=-oo 

On integrating by parts and using the periodicity conditions (39), we find that 


(42) 


J_ 

27T 



Z.^(u,v)e 


-inv 


dv = inS2j^(u) 






277 



Z^^(u, v)e~^*^^ dv = -n^f2j^(u) 


(43) 


We shall now show that each function can be determined as the solution of a cer- 
tain ordinary differential equation. To this end, multiply equations (27) and (33) by 
e“^ '^/2 t 7 and integrate both sides of the resvdting expressions with respect to v between 
-77 and 77. We find, on making the definitions 
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n 



V cosh H + cos H 
sinh H - sin H 





dv 


7 


n 



V sinh H - sin H 
cosh H + cos H 




> 


n = 0, ±1, ±2, . . . (44) 


and using equations (43), that must satisfy the ordinary differential equation 


Tn - 0, + 1, +2, . . . 

n” - {a + +n^ - 2in/3)f2j^ = < (45) 

Lo < u < Ug 


subject to the boundary conditions 


PoO^(O) + (s^ + inq^)S2j^(0) = 7jj(0) 

Ps"n<"s) + 


> n = 0,±l,±2, 


(46) 


where p , qi , s^, p„, q_, and s„ are defined by equations (34) and equations (21) de- 

O aJ O o o S 

fine a and /S as 


sin H sinh H 

(cos H + cosh H)(sin H - sinh H) du 

> 

^ ^ sinh H + sin H dH 

2 cosh H + cos H du J 


( 21 ) 


When h/d > 1/2, equations (44) to (21) (see eqs. (23) and (37) and ref. 3) can be re- 
placed by the approximate equation. 





7r 




g-inv 



n = 0,±1, ±2, . . . 


for ll > 1 (47) 

d 2 
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= 0 ,± 1 , ± 2 , 


n” - (n^ - 2in/3)J2j^ = ^ 


0 :< U < U 


for h>l 
d 2 


(48) 


J2;(0) + in[H(0) - Ipj^(O) = y^(0) 
n;(Ug) + in[H(Ug) - l]n„(Ug) = 


n = 0, ±1, ±2, . 


for ->- 
d 2 


(49) 


^ = 1^ forH>i 


2 du 


d 2 


(50) 


Finally, equations (25) and (42) show that, once the systems (45) and (46) have been 
solved, the surface displacement can be foimd from the relation 


V cosh H + cos H 
sinh H - sin H 



Oj^(u)e 


inv 


(51) 


or when h/d > 1/2 approximated by 

00 

” d 2 

n=-<=o 


(52) 


Thus, an exact solution for the sea-level elevation (and, since all other physical 
quantities can be expressed in terms of this, an exact solution to the complete flow prob- 
lem) can be obtained once the ordinary differential equation (45) with the boundary condi- 
tion (46) has been solved. This solution is given parametrically in terms of the variables 
u and V. These parametric variables are determined in terms of the physical variables 
X and y by the mapping 


z — w 


with z = X + iy and w = u + iv, which maps the lake or sea conformally into the rectan- 
gular region in the w-plane shown in figure 4 . The coefficients of the ordinary differen- 
tial equation depend on the particidar choice of the function H(u) . Thus it is impossible 
to proceed further in the general case. 

The procedures involved in obtaining a complete solution are best illustrated by con- 
sidering a particular case. In order to simplify the calculations it will be assumed that 
h/d > 1/2 and the approximate equations (47) to (50) will be used, A particular depth 
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distribution function H(u) will first be chosen. The ordinary differential equation will 
then be solved, and the surface displacement will be found as a function of u and v. A 
particular shape of the body of water and its island will then be chosen. Once this is 
done the relation between the physical variables x and y and the parametric variables 
u and V will be found. This will give the surface displacement as a function of x 
and y. 


SOLUTION OF ORDINARY DIFFERENTIAL EQUATION FOR 
PARTICULAR BOHOM TOPOGRAPHY 


We shall consider the case where the function H of u is given by 

H(u) = H+-iu 0<u<u^ (53) 

^s 

where is the constant depth at the mainland shore and + 6 is the constant depth 
at the island shore. On inserting equation (53) into equations (48) to (50) we obtain 


- [n^ - ~|f2 = r 

^ \ u J ^ " 


rn = 0,±1, ±2, . 


L.0 < u :< u^ 


(54) 


n^(0) + in(H„ - l)n„(0) = yjo) 


n' 


n' 


n = 0, ±1, ±2, . . . 


(55) 


fi;(Us) + in(H^ + 6 - l)fi„(Ug) = 


Before obtaining the solution to this boundary value problem, we shall first show that 


r(j(u)du = ro(Ug) - yq ( 0 ) 


(56) 
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To this end notice that equations (47) imply 



Hence it follows from equation (24) and the Cauchy-Reimann equations that 



2ir 


J_ 

2tt 





dz 

dw 


2 

du dv 





dv 



a(2)| du dv 

y J 9(u,v) 


On recalling that the rectangular region in the w-plane over which this integration 
is performed transforms into the interior of the body of water in the z -plane, it follows 
from the theory of transformation of multiple integrals that 



rQ(u)du = -i 
" 2tt 




dx dy 


(58) 
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where the integration is carried out over the entire siirface occupied by the body of 
water. Let C denote the distance measured along the boundary of this region; let 
S denote the distance measure along the boundary y of this region; and let the positive 
directions of C and S be counterclockwise and clockwise, respectively. Then applying 
the divergence theorem to equation (58) yields 




dS 


Using equation (32) gives 




Now 


and 


ia?|dc+-i(/ /?« L/'SiflWds (59) 

ldz| 2n f L \dw/ J |dz I 




Ids I =V(dx)2 + (dy)2 


dw 

dz 


\\dx) \dx) 


Hence 


dw 

dz 


)dS) = 




2 \ 2 
dx) 


\ 


.9x / \9x / \0x / 


On usii^ the Cauchy- Hi emann equations and noting that = “^y^x 


^1 |dS ( = y (u^ + "y + (Vjj dx + Vy dy)^ = ^(du)^ + (dv)^ 




But u is constant and equal to u on if. Hence, 
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|dS I = |dv ( 


dw 
dz 

or, taking into account the directions of transversal of the contours indicated in figures 4 
and 5, we see that 


dw 

dz 


dS = dv 


It also follows in the same way that 


dw 

dz 


dC = -dv 


Equation (59) can therefore be written as 



Finally, comparing this with equation (57) shows that equation (56) holds. 

Having established this result, we can now proceed to solve the boundary value prob- 
lem (eqs. (54) and (55)) for n = 0. In this case equations (54) and (55) become 


and 


0^(0) = ro(o) ' 


(60) 


(61) 


Integrating equation (60) gives 




FQ(t)dt + 
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where Z j is a constant of integration. On substituting this into the first boundary condi- 
tion (61), we find that Z j = ygCO). Hence, 

Oq- rQ(t)dt + r(j(0) (62) 

Equations (56) now show that the second boundary condition (61) is automatically satisfied. 
Notice that it would be impossible to satisfy both boundary conditions (61) if equation (56) 
did not hold. 

On the other hand, integrating equation (62) gives 

= r„(t)dt + y„(0)u 

where Zq is a constant of integration. An integration by parts yields 

= ^0 ^ n bx- t)EQ(t)dt (63) 

Since both boundary conditions (61) are satisfied, the solution (63) is indeterminate 
to within an arbitrary constant Zq. But this is as it should be since the n = 0 term in 
equation (52) is independent of v and as can be seen from equations (29) and (36) the 
original boundary value problem for ? involves only the partial derivatives of ^ and, 
hence, determines ^ only to within an arbitrary constant. 

We now return to the boundary value problems (54) and (55) and suppose that n ^ 0. 

Put 



where for definiteness the branch cut has been chosen to lie along the imaginary n-axis 
between 0 and i6/u . Then, the general solution to equation (54) is 

XU Xu 

J2j^ = Z^”^e’^ + ^ 2 ^^® ^ + -L / sinh[Xj^(u - t)]rj^(t)dt n = ±1, ±2, . . , (65) 

Jo 
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where and are constants of integration. Substituting this expression into 
boundary conditions (55) shows that and satisfy the equations 

+ ln(H„ - I)] - - in(H„ - 1)] > rJO) 


+ ln(H^ + e - ' “<®o + 6 - l)]e 


- >'n<“s> - 


r"= 

if 

n Jq 


cosh X^(Ug - t) + in(HQ + 6 - l)sinh Xj^(Ug - t)]rj^(t)dt 


Solving these equations for and Z^*^^ gives 


zl^^K 

= 

Xn - in(Ho - D 


r(n)j 

2 ^ 


n 


\ + ‘"(H„ - X) 


'Jn- 


rn(»)Kn 






where we have put 

K„ . . ln(H^ - + ln(H„ + 6 - 


X„u 

n s 


[X„ + in(H^ - l)][Xj^ - in(HQ + 5 _ i)]e ® n = ±1, ±2, 


(66) 


Jn = ’^n^^s) - >'n( 0 ) 


Xn - + 6 - 1 ) I ^ 


Xn - in(H^ - 1 ) 


— I w„(t-u^;H 
2x I n so 

n Jq 


+ 6)rj^(t)dt n-±l,±2, 


(67) 


21 



and 


w. 


0 ^ [A„ - in(? - + [X^ + in(? - l)]e n = ±1, ±2, . . . (68) 


Substitutii^ these results into equation (65) now shows that 




-X„u 

yn(0)e ^ 


Xn - in(H^ - 1) 


if 


+ — I sinh[Xj^(u - t) ]r^^(t)dt 

'O 


n = ±1,±2, . 


(69) 


Substituting equations (63) and (69) into equation (52) and noting that [fij^(u)j* = 
gives 



^r. 

Jn 




sinh[X^(u - t)]rj^(t)dt 


,inv 


(70) 


This gives the general solution to the problem in terms of the integrals Tj^(u) and yj^(u), 
which are defined in equations (47). It can be seen from equation (17) that, in order to 
evaluate these integrals, the complex wind stress and the function z — w, which 
maps the interior of the body of water into the rectangular region in the w-plane (shown 
in fig. 4), must be specified. For lakes a good approximation is to take the wind stress 
as a constant. We shall therefore work out specific results only for this case and in the 
remainder of the report we wiU assume that is a complex constant. The function w 
can be found by the techniques of conformal mapping once the shape of the lake and island 
are specified. 
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Reduction of Integrals and In the Case of Constant Wind Stress 

When the wind stress is equal to a complex constant, it is possible to reduce the 
integrals (eq (44)) for and to a form in which the only integrals that appear are 
contour integrals which depend only on the function w and not on the depth distribution 
H. In order to simplify the calculations we shall restrict ourselves to the approximate 
form (eq. (47)) with a given by equations (17) and (19) for h > id. 

It is convenient to introduce the functions z and z* as independent variables. 

Then, since w is analytic, w is a function of z only and not of z* , and w* is a func- 
tion of z* only and not of z. In addition, since 

z = X + iy 


and 


z* = X - iy 


Then 


Hence, 



da _ 

1 

'aa(l) , 0a(2) 

i 

+ 

'aa(2) 


dz 

2 

dx dy 

2 

8x 

9y 


az 2 L* y J 


Since z is a function of w only and not of w this can be written as 


^ y Vaw dz ) 


(71) 
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but, when is a constant, it follows from equations (17) and (19) that a is 
of u only and independent of v. Hence, 


and 


9g _ 

1 ^ 

3w 

2 du 

Tda dwl 

= i[(' 

(_du dz J 

2LV( 


Also, 


+aY-)l 

\dw/ 2 L Vdw/ Vdw/J 


After inserting eqviations (72) and (73) into equation (47) and recalling that 
independent of v and that 


dz 

dw 


dw \dw/ 


we get 


Tn(u) = ± 
^ 2 


J. f ^ dv f 

\du/ 2tt J ^ dw \du/ 2u J_^ \dw, 

»=i L*J- /" e-*”' i dv + a ^ 

J-1 J-t 


* 

^ dv 


or, since w = u + iv. 


function 


(72) 


(73) 

CT is 
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i 

I 

it 






n = 0,±1, . . . 


where we have put 


- > i (74) 
d 2 


I 


n 


2ir 



e-nw^dv = JL 

dw 27Ti 



e-”^^dw 

dw 


It can be seen from figures 3 and 4 that by making the change of variable w = In T we 
transform this integral into an integral about a closed contour, C^, for example, in the 
T -plane. Thus 


I 


n 


1 

27ri 



J- ^ dT 
rr.n dT 


(75) 


where, since dz/dT has no singularities inside the annular region between the circles 
of radii and 1, the contour can be taken as any closed contour lying within this 
region and encircling the origin. 

Now it is well known from the theory of functions that dz/dT must have a Laurent 
series expansion about the origin in the T-plane that converges in the annular region. It 
is easy to show that Ij^ is simply the coefficient of T^”^ in this series. In the general 
case an identical analysis beginning with equation (47) shows that 


/cosh H + cos H 

7d 


aV-nu j* ■ 

(- 

sinh H - sin H 

\d 

u/ ” Vd 

u/ "" 


yu) = i £ 

” 2 V sii 

V (u) .1 >lnhH-Bln ;i r 
2 Y cosh H + cos H L 


* „nu T ^-nu X* 
e Iji ^ 




n = 0,±1,±2, 


It follows from figures 1 and 2 and equation (75) that 



dz = 0 
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Hence, 


ro(u) = = 0 

It follows from equations (17) and (19) that 

c = l 2 e-[( W)/2]H A for £ > i 

gd \ J d 2 

Hence, 

^ = _ 27riT^ (1 ^ j)g-[(l+i)/2]H m for ]L>1. 
du gd du d 2 


(76) 


(77) 


(78) 


Simplification of Solution for Surface Displacement 
in the Case of Constant Wind Stress 


We shall now use the results obtained in the previous section to simplify solution (70) 
for the case where the wind stress is constant. 

To this end notice that equations (63) and (76) show that 


S2q = Zq = constant 


Now for n 0 it follows from equation (64) that 


* 

r. 

n -n 


and therefore equations (66) and (68) show that 


♦ 

n -n 


= -K „ 
n -n 


(79) 


(80) 


(81) 

(82) 
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Equations (67), (69), and (74) show that we can now write 


J =ll* 
n 2 ^ 


+ -- In 
2 n n 


( 2 ) 




( 1 ) 

n 



(83) 


where 




Xn - in(H^ -H 6 - 


Xn - in(H^ - 1) 


n c/n 


+ 6) ^ e-"^ dt 
dt 


(84) 


J® = /(u^)a”"=-a*(0) 


X„ - ln(H^ + e - l)l^-Vs 


- 1 ) 


'iT / “n«-V«o 

n Jq 

(1)_4^\ / g(0)e , 1 f , 

" K ” ’ ° X - in(H - 1) X / sinh[X (u - t)]^e‘ 
n n ' o ' n dt 


+ 6)/^') e”* dt 

Vdt 


(85) 


and 


t(2) 

(2) ^ jj X _ CT J 0 ) e ^ / sinh[X (u -t)](^] 

K„ " ’ o X„-in(H -1) X„ / AdtJ 

n n o ' n 


e”*dt 
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Equations (66), (80), and (82) now show that 


o(0)K 


-n 




^ [^-n - - 1) ^[^-n + ~ ] 


therefore, equations (66), (68), and (80) to (82) show that 


_ M 3fe 

l^nj n = ±l,±2. 


(86) 


Recalling that (see remarks preceding eq. (70)) ~ ^ n using equation (79) we find 


? = Zq + 2 /?e ^ a^(n)e 


inv 


n=l 


Substituting equation (83) into this expression gives 


n=l n=l 

Since for any complex number M, Re M* = ReM, this becomes 




z 

n=l 




-inv 


Finally, using equation (86) and the fact that Iq = 0, shows that 


^-Iq+Rb^ 


n=-oo 


(87) 


After integrating by parts equation (84) becomes 


4‘> 


(t(O)K. 


n 


2Xj\n-i"^«o- 2X. 
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Substituting this into equation (85) and performing another integration by parts yield 




( 1 ) _ 


n 


n n 




dt 



Substituting equation (53) into equation (77) and using the result in equation (88) , gives 
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Since 
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and 




sinh Xj^(u - t)]dt 



It follows that 


„a) , ^ jj(o) 

n gd " 


(89) 


where we have put 

^ - "n<« - “s>Ho + «“n(Ho>] 
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and 


Since 


- 


n5(l + i)r2 + in 




- 1 



.2 _ (1 + i)6 



substituting equation (89) into equation (87) gives 




gd 


? I «' 


(u)e 


inv 


n=-oo 


(90) 


where 




nu. 




|nu+(l!:l]H(u) 


2e 


2n + (1 + i) (n + 1 A] 


(91) 


2 u_ 


Wn(u; 0 = [\ - in(? - l)]e " + [\^ + in(^ - l)]e ” 


( 68 ) 


K„ - [X^ - in(H^ - 1)][X^ + in(H^ + 6 - l)]e 


Xu 
n s 


- [^n + i"(«o - + 5 - l)]e 


( 66 ) 
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2n(2 + ip 


(92) 


Mn(?) 


2n + ( 1 + i) 






- 1 


I =J_ (f J-^dT (75) 

27Ti Jq dT 

This gives the solution to the problem for the case of a constant wind stress. Numerical 
results can be obtained once the constants 1^^ are specified. These constants depend 
only on the shape of the lake and island, and it is only through these constants that the 
specific lake and island geometry enter the problem. They can be determined once the 
mapping z — T of the region occupied by the body of water into the annular region of 
figure 2 is known. The constants 1^^ are then determined either by evaluating the con- 
tour integral (75) or by finding the coefficients of the Laurent series expansion of the 
mapping about the origin in the T -plane. 

Since w = In T, a knowledge of the mapping z — T is equivalent to a knowledge of 
the mapping z — w, which transforms the cut region occupied by the body of water shown 
in figure 5 into the rectangular region in the w-plane shown in figure 4. This mapping 
can be used to express the variables u and v in the w plane as functions of x and y. 
Since equation (90) gives the surface displacement as a function of u and v, a knowledge 
of the mapping z — w allows us to express the surface displacement parametrically as a 
function of position in the physical plane. In order to illustrate these ideas a complete 
numerical solution will be worked out for a specific geometry. 


Solution of Complete Boundary Value Problem for a Circular Lake Containing an 
Eccentricity Located Circular Island With an Imposed Constant Wind Stress 

In order to find the mapping z — w the specific shape of the body of water and its is- 
land must be specified. We shall, therefore, suppose that the body of water is a circle 
with radius r and that the island is also a circle whose radius is rg. However, the 
center of the island and the center of the lake do not necessarily coincide. The plan view 
of the lake and island are shown in figure 6. As disciissed, the slit is included in the fig- 
ure, and with no loss in generality we have chosen the x-axis to lie along the line of cen- 
ter of the island and the lake. In order to find the analytic function w with the appro- 
priate properties it is necessary to map the cut eccentric annular region shown in 
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figure 6 into the rectangular r^ion in figure 4 in the manner indicated in the figures. 
This is best accomplished by first finding the mapping z — T, which transforms the ec- 
centric annular region in the z -plane (fig. 6) into the concentric annular region in the T- 
plane as shown in figure 3. But this mapping is given in reference 4 (p. 287) as 

T = — ( 93 ; 
az - r 

where a is related to the coordinates and x^^ of the island shown in figvtre 6 by 



with a > 1 and where the radius of the larger circle in the T -plane (shown in figs. 2 
and 3) must be equal to 



Since w = In T it follows from equation (93) that 


w = u + iv = In 



z = X + iy 


(95) 


(96) 


Hence, for each value of the point x,y in the region occupied by the body of water in 
figure 6 we can compute a point (u,v) in the rectangular region of figure 4. For this 
value of u and v, we can find a value of the surface displacement ^ from equation (90). 
This determines ^ as a function of position in the physical (z -plane) plane. 

It only remains to determine the constants Ij^. To this end we solve equation (93) 
for z to obtain 


z 


= r 


a - T 
1 - aT 
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Differentiating with respect to T gives 


^ = r ~ ^ 
(1 - aT)^ 


On inserting this into equation (75) we obtain 


I 


n 


r(a2 - 1) J_ 

2vi 



i dT 

t "(1 - aT )2 


n = ±1, ±2, . . . 


(97) 


where is any contour encircling the origin and lying between the unit circle and the 
circle of radius in the T -plane (fig. 2). This integral can easily be evalmted by the 
method of residues. Since a > 1 the pole of order two of the integrand at the point 
T = 1/a lies within the unit circle and therefore within the contour C^. In addition for 
n > 0, Ijj has a pole of order n within at T = 0, but this pole will not occur when 
n < 0. 

2 

In order to evaluate the residue at T = 0, we expand (1 - aT)“ in a Taylor series 
about T = 0 to obtain 


oo 

i =V (m + l)(aT) 


Hence the integrand has the Laurent expansion about T = 0: 


— ^ — =y 

T“(l-aT)2 ^ 


(m + l)a“T“"" 


1 n 1 

The coefficient of T~ is na"“ ; hence, at T = 0 the residue of the integral in 


equation (97), /?&f.(T = 0) is 


/?eii(T = 0) =< 


na 


n-1 


for n = 1, 2, . . . 
for n = 0, -1, -2, . . . 


35 


The residue Red- (T = i) at T == i is 
V a/ a 

= -V 

V ay \dT ^2^n/ 


= -na”“^ for n = 0, ±1,±2, 


T=l/a 


Hence, on taking the sum of the residues and using the residue theorem, we find that 




-r(a^ - l)na’^“^ for n= -1,-2, . . . 


for n = 0, 1, 2, . . 


=f: 


~ for n = 1, 2, . . . 

„n+l 


(98) 


for n = 0, -1, -2, 


Equation (98) gives the values of the constants that are needed for equation (90) to com- 
pute the surface displacement. 


Computation of Flow Field from Surface Displacement 

Having solved for the surface displacement, the other quantities of physical interest 
can be calculated. These quantities are the complex horizontal velocity U + iV which is 
given in terms of the svirface displacement by equation (10) and the depth averaged veloc- 
ity Q, which is given in terms of this surface displacement by equation (14). 

Since the surface displacement is computed first as a function of u and v, it is con- 
venient to express the other quantities as functions of u and v and then use equa- 
tion (96) to relate them to x and y. 

In order to accomplish this, we use equation (71) to write 


0Z 


but, since w is anal 5 d;ic, w is a function of z only. Hence, 
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Inserting this together with equation (19) into equation (10) now yields 



( 100 ) 


Upon noting that T = e^, we find that 


dz _ dz dT _ ^w dz 
dw dT dw dT 

It now follows from equation (97) that the term (dw/dz) in equation (100) can be ex- 
pressed as a function of u and v by 



( 101 ) 


Thus, eqiiations (90), (100), and (101) give u + iv as a function of u and v and 

In a similar way we can substitute equations (19) and (99) into equation (14) to ex- 
press Q as a function of u and v: 


= a - |l + i[l - H(u)]j(?^ + (102) 

However, instead of working with the volume flow itself, it is more convenient for plot- 
ting purposes to find the stream function. In reference 1 it is shown that Qj and Qg 
satisfy the continuity equation 


9Qi 9Qp 
— ^ + — ^ = 0 
dx By 


Hence, there exists a stream function if/ such that 


ay 




Qt2 = - — 
Bx J 


(103) 


Hence 


Q = ^ +Ujy 

By Bx i ^ ^ i i \dz / i \9u 9vAdz / 

Inserting this into equation (102) dividing through by 1/i (dw/dz) and taking the real 
part gives 


27rf 


c Bx}/ _ 

gd 3u 



(104) 


Since the normal volume flow vanishes at the boimdaries the boundaries themselves must 
be streamlines. Hence xj/ must be constant along each boundary. Since x}/ is only de- 
termined by equation (103) to within an arbitrary constant, we can set it equal to zero 
along the mainland coast 12 in figure 6. Since this line maps into the line u = 0 in the 
w-plane (fig. 4) it follows that 


X}/(0,V) =0 -7T < V < 7T 

Hence, equation (104) can be integrated along a line of constant v to obtain 


— -x}/{u,v)= /?e i / a( — ^ du + / 

gd Jq Vdw/ 

which together with equations (17), (90), and 
of u and v. 


[1 - h(u)]lidu + -L 
0u 9v 



(101) allow us to compute 


^}/ 


du (105) 

as a function 
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DISCUSSION OF RESULTS 


A method for obtaining solutions to a certain class of partial differential equations 
in a doubly connected domain is given. This class of equation includes the eqviation for 
the wind -driven circtilation in shallow seas and lakes for a large class of bottom toixag- 
raphies. The technique is used to reduce the problem of finding solutions to the shallow 
sea equation for a closed body of water with an island to the problem of solving an ordi- 
nary differential equation. An approximate form of this ordinary differential equation, 
which is valid when the depth of the body of water is greater than one half of the Ekman 
depth, is solved for a particular bottom topography and a constant wind shear stress. 

The results presented herein are for a circular lake containing a concentrically or 
eccentrically located circular island with a logarithmic bottom topography. Two classes 
of bottom topography are included: One where the lake is deepest at the mainland shore 
and the other where the lake is deepest at the island shore. The assumption that the 
water depth is greater than one half the Ekman thickness is quite realistic. For instance 
in Lake Erie, which is the most shallow of the Great Lakes, only some 15 percent of the 
total bottom in the central and eastern basin is less than 1/2 the Ekman thickness for a 
10 meter per second (22 mph) wind. 

The surface displacement of the body of water is given in terms of the intermediate 
variables u and v by equation (90) with the aid of equations (91), (68), (66), (92), and 
(98). The horizontal velocity as a function of u and v and depth ^ is given by equa- 
tion (100) and the stream function i// as a function of just u and v is given by equa- 
tion (105). Since the physical coordinates x and y of the horizontal plane of the body 
of water are related to the intermediate variables u and v parametrically by equa- 
tion (96), the surface displacement, the horizontal velocity, and the stream function are 
known as functions of x and y. The solution for U, V, and ij/ depend on the wind 
shear stress t^, the friction depth d, the bottom depth h(x,y), the basin radius, the 
island diameter (x„ - Xv.)/r, and the eccentricity of the island defined as (x„ + x. )/2r. 

ci U I W I 3. D 

The magnitude of the wind stress | t | enters only as a normalizing parameter in the 
solutions. Table I contains a list of the various solutions for which results are given 
herein. 

The bottom topography with an island diameter of 0. 5, island eccentricity of 0, 
h/d = 1. 0 at the mainland shore, and h/d = 0. 5 at the island is shown in figure 7(a) . 

This topography includes the main effects due to an island. In an actual basin without an 
island, the depth is smallest at the mainland shore and largest at some central interior 
point. Islands create local areas in the interior of the basin where the depth becomes 
shallow. The topography shown in figure 7(a) will therefore give the local effect due to 
the island. 
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TABLE I. - CIRCULAR BASIN SOLUTIONS 


Figure 

number 

Depth ratio, h/d 

Is 

land 

tricity, 

\)A2r) 

Island 
diameter, 
(Xa - x^,)/r 

ryn 

At mainland 
shore 

At island 
shore 

eccen 

(Xa + 

7 

1.0 

0. 5 

0 

0. 50 

i 

8 

1.0 

.5 


.5 


50 

i 

9 

1.0 

.5 


.5 


50 

1.0 

10 

1.0 

. 5 

0 


25 

i 

11 

3.0 

. 5 

0 



50 

i 

12 

. 5 

1.0 

0 



50 

i 

13 

.5 

3.0 

0 



50 

i 

14 

.5 

3.0 

0 



25 

i 

15 

. 5 

1.0 


.5 


50 

1.0 

16 

. 5 

1.0 





i 

17 

. 5 

3.0 





1.0 

18 

.5 

3.0 





i 

19 

3.0 

.5 





1.0 

20 

3.0 

. 5 
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The surface displacement caused by a wind in the positive y-direction for the bottom 
topography depicted in figure 7(a) is shown in figure 7(b). Gedney (ref. 2) has shown that 
the solution for a constant depth basin is the same whether or not an island is present. 

As is well known, the surface displacement solution for a constant depth basin is a plane 
inclined to the horizontal. The deviation of the surface displacement shown in figure 7(b) 
is then the effect of the variations in bottom depth. Although the deviations of the surface 
displacement from an inclined plane are small, they will have major effects in the local 
velocities. 

The horizontal volume flow stream function is shown in figure 7(c) for the case under 
discussion. The stream function plot consists of two gyres; the gyre to the right of the 
wind is rotating clockwise and the one to the left of the wind is rotating counterclockwise. 
The two gyres are separated by a dividing streamline, which in this case has a value of 
zero. The value of x(/ on both the island and mainland shores is zero, and the zero 
streamline in the interior has been twisted relative to the wind in a clockwise direction. 

As shown in the ij/ plots in reference 2, the dividing streamline for a dish shaped 
basin whose depth increases monotonically from the shore to some interior point always 
runs through the maximum depth point. If the depth decreases monotonicaUy from the 
shore to an interior point, the dividing streamline would run through the minimum depth 
point. In the case of the island shown in figure 7(c), the zero streamline runs through 




the minimum depth point interior to the mainland shore which happens to be the island 
boundary. If the island was not placed at the minimum depth point, the value on the 
island boundary would assume some other value; that is, the path of the streamline divid- 
ing the two gyres is determined primarily by the extremes in the bottom topography. 

The horizontal velocities at ^/d = 0, -0. 125, -0. 250, -0. 375, and -0. 500 are shown 
in figures 7(d) to (h). The magnitude of each velocity vector plotted can be determined 
from the scale included on the plot. The origin of the velocity vector is the position at 
which velocity is actually occurring. The velocities at the surface are skewed to the 
right of the wind mainly because of the Coriolis force. Local perturbations in the sur- 
face velocities due to the volume flow gyres are evident. The acceleration of the veloc- 
ities arotmd the island and the deceleration near x = -1 and +1 can be directly attributed 
to the volume flow gyres shown in figure 7(c). At ^/d = -0. 125 mass is still being trans- 
ported in the direction of the wind but an increasing amoimt is being transported to the 
right of the wind. The flow acceleration around the island is still evident. At 
^/d = -0. 25 the flow pattern is very similar to that of the integrated volume flow shown 
in figure 7(c). Return flow opposite to the wind direction is now occiorring. The flow 
near the island is predominately tangent to it. At ^/d = -0. 375 and -0. 50 there is much 
return flow opposite in direction to the top surface layer flow. The flow below 
^/d = -0. 50 is very similar in direction but decreases in magnitude from that shown at 
^/d = -0. 50. This general flow pattern where mass in the top siirface layer is trans- 
ported in the direction of and to the right of the wind and returned in the opposite direc- 
tion in the bottom layer is the dominant pattern in shallow lakes. The variations in this 
general pattern are due to the particular bottom topography. For this particular case 
where the bottom slopes upwards at the island an acceleration of flow occurs aroimd the 
island. 

The effect of offsetting the island relative to the center of the basin is shown in fig- 
tire 8. The bottom topography, as shown in figure 8(a), has the same values at the main- 
land and island as the previous case but the center of the island is at x/r = 0. 5. The 
wind direction and island diameter are the same for the two cases shown in figures 7 
and 8. The effect of the eccentricity on the stream function is clearly shown in fig- 
ure 8(c). The zero streamline still runs around the island shore since this is still the 
shallowest point in the lake. The gyres though no longer symmetric stiU have the same 
overall behavior as the concentric island case. Comparing the velocity plots shown in 
figures 8(d) to (h) with the ones in figure 7 shows that the flow patterns are essentially 
the same type for both cases. Making the island eccentric has the main effect of locally 
contracting and expanding the patterns for the concentric island case. 

Figure 9 gives the results for an island diameter of 0. 5, an eccentricity of 0. 5, 
h/d = 1. 0 at the mainland shore, h/d = 0. 5 at the island shore, and with the wind along 
the x-axis. In the concentric island case the bottom topography is independent of the 
angle about the origin of the x,y axes. As a result, the solutions for any two wind 
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directions are the same except rotated by the difference in wind angles. Therefore, it is 
possible to compare the results shown in figure 7 with those in figure 9 to determine the 
effect of an eccentric island when the wind is along the x-axis. As in the case where the 
wind is in the direction of the y-axis, the effect of locating the island eccentrically when 
the wind is along the x-axis is to cause a local contraction and expansion of the patterns 
which occur in the concentric island case. 

The effect of change in island size can be seen by comparing figure 7 and figure 10. 
The results shown in these figures are for the same parameters except the island diam- 
eter is 0. 5 in the case of figure 7 and 0. 25 in the case of figure 10. The bottom topog- 
graphies for these two cases are very nearly the same when plotted against distance from 
the island shore. As a result, for points equidistant from the island shore the velocities 
in the two cases are very nearly the same in magnitude and direction. This just illus- 
trates again the importance of bottom topography. 

The combined effect of making the basin deeper and the bottom slopes greater is 
demonstrated in figure 11. Here h/d = 3. 0 at the mainland shore, and h/d = 0. 50 at the 
island as shown in figure 11(a). Figure 11(b) shows first that the average surface dis- 
placement slope is much smaller than for the case shown in figure 7(b). This, of course, 
is because the basin shown in figure 11 has a greater average depth than the basin shown 
in figure 7. However, the surface displacement pertvirbation due to the change in bottom 
depth are much greater in figure 11(b) than figure 7(b). This effect is due to the larger 
bottom slopes in the case shown in figure 11. 

The effects of larger depths and bottom slopes on the stream fimction can be deter- 
mined by comparing figures 11(c) and 7(c). First, the zero streamline is twisted clock- 
wise relative to the wind by a much greater angle in the case shown in figure 11(c). Sec- 
ond, the value of the stream function gradients are much greater in figure 11(c), creating 
larger volume flows. 

The stream function effects just mentioned produce the major differences between 
the flow velocities in figures 11(d) to (h) from those in figures 7(d) to (h). Because the 
zero streamline in figure 11 is more perpendicular than parallel to the wind, the flow 
near the surface and bottom are more alined to the horizontal axis than in figure 7. Also 
in figure 11 there is a greater acceleration of velocities arovmd the island. 

When the sign of the bottom slope (0h/9x, 9h/3y) is reversed, we can expect a signif- 
icant change. In figure 12, the results are given for h/d = 0. 5 at the mainland shore, 
h/d = 1. 0 at the island, eccentricity = 0 and island diameter = 0. 5. The average svtrface 
displacement slope in figure 12(b) is greater than in figttre 7(b) because the average depth 
in the case shown in figitre 12 is less than in the case shown in figure 7. There are dif- 
ferences in the deviations from a constant slope surface displacement between the two 
cases although they are not easily discernible. 

The effect of changing the sign of the bottom slope can be more readily seen by com- 
paring the stream function plots in figures 12(c) and 7(c). When for x < 0 9h/9x is 
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positive as in the case of figure 12 the zero streamline is rotated counterclockwise in- 
stead of clockwise as in figure 7(c). The resulting effect on the flow velocities are easily 
seen in figures 12(d) to (h). The rotation of the zero streamline relative to the wind in a 
counterclockwise direction is characteristic of many basins where the depth increases 
from the mainland shore to an interior point. For example reference 2 shows this effect 
for a rectangular basin without an island. As we found previously, the amount of rotation 
of the zero streamline depends on the depth of the basin and bottom topography slopes. 

The streamline patterns for h/d = 0. 5 at the mainland and h/d = 3. 0 at the island are 
shown in figures 13(c) and 14(c), respectively, for island diameters of 0. 5 and 0. 25. 

Near the islands, the zero streamline for these cases have been twisted counterclockwise 
considerably more than in the case shown in figure 12(c). This local increase in rotation 
is due to the large bottom slopes near the islands as shown in figvtres 13(a) and 14(a) as 
compared with the ones in figure 12(a) . 

Other cases noted in table I are also included for the interested reader but are not 
discussed. 


CONCLUDING REMARKS 

An analytical solution has been obtained for the wind -driven circulation in a shallow 
lake containing an island. The effect of an eccentric circular island in a circular lake 
where the bottom depth is decreasing from the mainland shore to the island is shown. 
This island effect is contrasted with the more usual case where the bottom depth in- 
creases from the mainland shore to some interior point. The island effect is shown to 
produce a completely different volume flow pattern (here the volume flow refers to the 
vertically integrated velocities). These volume flow patterns result in unusual velocity 
patterns that include the acceleration of flow near the island shore. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 6, 1971, 

129-01. 
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APPENDIX - SYMBOLS 


A 

a 

B 

C 

c. 

D 

d 

E 

F 


f 

f 


c 


H 

H 

1 

h 



ja),j(2) 


n 


111 I 

Mn(?) 


n 

A 

n 


" l >^2 

Po’^-Ps’^ls 


function of h/d given by eq. (2) 
mapping parameter defined by eq. (94) 
function of h/d given by eq. (2) 

function of h/d defined in eq. (5) or distance along mainland shore 

mainland shore line 

contour of integration in T -plane 

function of h/d defined in eq. (5) 

Ekman friction depth, 7r(2iV'f^)^^^ 
functions of h/d defined in eq. (6) 
functions of h/d defined in eq. ( 13) 
defined by eq. (28) 
function defined by eq. (22) 

Coriolis parameter 
acceleration due to gravity 
defined in eq. ( 19) 

value of H at lake or sea mainland shore 

depth of lake or sea 

defined in eq. (75) 

imaginary part 

defined in eq. (67) 

defined in eq. (84) 

defined by eq. (66) 

constants of integration 

constant defined in eq. (92) 

integer 

unit normal to boundary of sea or lake 
X- and y -components of n 
constants defined in eq. (34) 
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Qi,Q2 


s , s 
o’ s 


rn(u),ynN 



dummy variable of integration 

x-component of horizontal velocity 

harmonic function of x and y 

u-coordinate of boundary of rectangle in w-plane 

y-component of horizontal velocity 

harmonic conjugate of u 

complex function u + iv 

coordinate of lake surface plane 

x-coordinates of Island boundary points 

coordinate of lake sxirface plane 

transformed dependent variable 

X + iy 

defined in eq. (21) 

Fourier coefficients defined in eq. (44) 
parameter in depth distribution, eq. (53) 
displacement of sea or lake surface 
dummy variable 
eigenvalue defined in eq. (64) 
coefficient of vertical eddy diffusivity 
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w 


w w 

n ’ n ’ n 
w^(u; C) 


Subscripts: 


vertical distance measured upwards from water surface 

. ia(2) 

defined in terms of wind stress by eq. (9) 

. lr» 

X, y components of wind stress divided by density 
stream function for the volume flow Q 
Fourier coefficient defined by eq. (41) 
defined by eqs. (85) and (91) 
defined in eq. (68) 


u,v partial derivative 

X, y partial derivative 


Superscript: 

* 


denotes complex conjij^ate 
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(b) Surface displacement. 


152 

)5 


d shore; h/d = 0. 5 at island; T’'y|T*''( = i. 











(d) Horizontal velocities at f/d = -0. 125. 



(e) Horizontal velocities at 
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(f) Horizontal veixitles at £/d = -0. 250. (g) Horizontal veixities at f/d = -0. 275. 
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(h) Horizontal velocities at f/d = -0. 500. 


Figure 9. - Concluded. 
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(e) Horizontal velocities at f /d = -0. 125. 
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(f) Horizontal velocities at f /d = -0. 250. 
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(g) Horizontal velocities at f/d = -0.375. 
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(h) Horizontal velocities at f /d = -0. 500. 
Figure 12. - Concluded. 




(a( Bottom topography. 


(b) Surface displacement. 



(c) Nondimensional stream function 

Figure 14. - Circuiar basin circuiation. Isiand eccentricity = 0; isiand diameter = 0. 25; h/d = 0.5 at main iand shore; hfd = 3.0 at island; 
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(f) Horizontal velocities at f /d = -0. 250. (g) Horizontal velocities at |/d = -0. 375. 



Figure 18. - Concluded. 
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(al Bottom topography. 



(b) Surface displacement. 


(c) Nondimensional stream function 2?rf(.^'/|T''f|. 

Figure 19. - Circular basin circulation. Island eccentricity' 0.5; island diameter = 0. 5; h/d » 3.0 at mainland shore,- h/d = 0. 5 at island; 
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